Dynamical response of the "GGG" rotor to test the Equivalence Principle: theory, 
simulation and experiment. Part I: the normal modes 
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Recent theoretical work suggests that violation of the Equivalence Principle might be revealed 
in a measurement of the fractional differential acceleration r\ between two test bodies —of different 
composition, falling in the gravitational field of a source mass— if the measurement is made to the 
level of rj ~ 1CU 13 or better. This being within the reach of ground based experiments, gives them 
a new impetus. However, while slowly rotating torsion balances in ground laboratories are close to 
reaching this level, only an experiment performed in low orbit around the Earth is likely to provide 
a much better accuracy. 

We report on the progress made with the "Galileo Galilei on the Ground" (GGG) experiment, 
which aims to compete with torsion balances using an instrument design also capable of being 
converted into a much higher sensitivity space test. 

In the present and following paper (Part I and Part II), we demonstrate that the dynamical 
response of the GGG differential accelerometer set into supercritical rotation —in particular its 
normal modes (Part I) and rejection of common mode effects (Part II)— can be predicted by means of 
a simple but effective model that embodies all the relevant physics. Analytical solutions are obtained 
under special limits, which provide the theoretical understanding. A simulation environment is set 
up, obtaining quantitative agreement with the available experimental data on the frequencies of the 
normal modes, and on the whirling behavior. This is a needed and reliable tool for controlling and 
separating perturbative effects from the expected signal, as well as for planning the optimization of 
the apparatus. 

PACS numbers: 04.80.Cc, 07.10.-h, 06.30.Bp, 07.87,+v 



I. INTRODUCTION 

Experimental tests of the Equivalence Principle (EP) 
are of seminal relevance as probes of General Relativity. 
The Equivalence Principle is tested by observing its con- 
sequence, namely the Universality of Free Fall, whereby 
in a gravitational field all bodies fall with the same accel- 
eration regardless of their mass and composition. They 
therefore require two test masses of different composi- 
tion, falling in the field of another 'source' mass, and a 
read-out system to detect their motions relative to one 
another. An EP violation would result in a differential 
displacement of the masses in the direction of the source 
mass, which cannot be explained on the basis of known, 
classical phenomena {e.g. tidal effects). 

The landmark experiment by Eotvos and collabora- 
tors Q has established that a torsion balance is most 
well suited for ground tests of the EP, thanks to its inher- 
ently differential nature. With the test masses suspended 
on a torsion balance they improved previous pendulum 
experiments by almost 4 orders of magnitude, showing 
no violation for r\ larger than a few 10~ 9 0. Several 
decades later, by exploiting the 24-hr modulation of the 
signal in the gravitational field of the Sun, torsion balance 
tests have improved to 10~ n and then to 10~ 12 
More recently, systematic and very careful tests carried 
out by Adelberger and co-workers using rotating torsion 



balances have provided even more firm evidence that no 
violation occurs to the level of 1CU 12 0,0. 

The relevant theoretical question for Equivalence Prin- 
ciple tests is: at which accuracy level a violation, if 
any, is to be expected? In an earlier work by Damour 
and Polyakov, based on string theory and the existence 
of the dilaton 6], rj values at which a violation might 
be observed have been determined to be in the range 
10 -18 < r\ < 1CU 13 . Fischbach and coworkers have 
derived a non-perturbative rigorous result, according to 
which a violation must occur at the level of rj ~ 10 -17 , 
due to the coupling between gravity and processes of v— D 
exchange which should differently affect masses with dif- 
ferent nuclei. More recent work suggests, in a new 
theoretical framework for the dilaton, that a violation 
might occur already at the level of 77 ~ 10~ 12 — 10 -13 , 
depending on the composition of the masses. 

While an r\ ~ 10~ 13 , and perhaps smaller, should be 
accessible with rotating torsion balance experiments on 
the ground, a sensitivity as high as rj ~ 10~ 17 could be 
achieved only by an experiment flying in low Earth or- 
bit, where the driving acceleration is up to three orders 
of magnitude larger. Specific instruments have been de- 
signed to carry out such an experiment in space: STEP, 
Microscope and "Galileo Galilei" (GG) [|-[l3- They 
share two features: that the test masses are concentric 
cylinders, and that rotation of the spacecraft provides 
signal modulation at frequencies higher than the orbital 
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one. 

GG is peculiar in that it spins around the symme- 
try axis and is sensitive to relative displacements in the 
plane perpendicular to it: cylindrical symmetry of the 
whole system and rotation around the symmetry axis al- 
low passive attitude stabilization of the spacecraft with 
no need of a motor after initial spin up to the nominal 
frequency (typically 2 Hz). The planar (instead of linear) 
sensitivity of the instrument is also a crucial feature for 
allowing us to rotate at supercritical speeds, i.e. faster 
than the natural frequencies of the system. Faster rota- 
tion means modulation of the signal at higher frequency 
and therefore a reduced 1/ / noise (for 1// noise see e.g. 
the website maintained by W. Li [jjj). GG differs from 
the other proposed space experiments also in that the 
test masses are suspended mechanically. We find that 
in absence of weight, as it is the case in space, mechani- 
cal suspensions too can provide extremely weak coupling, 
with the additional advantage to electrically ground the 
test masses. 

The GG design naturally allows us to build and test 
a full scale 1-g version of the apparatus: by suspend- 
ing the instrument on a rotating platform through its 
spin/symmetry axis, the sensitive plane lies in the hori- 
zontal plane of the laboratory where a component of an 
EP violation signal might be detected, similarly to a tor- 
sion balance experiment. "Galileo Galilei on the Ground" 
(GGG) 0, 0j is primarily a prototype for testing the 
main novel features of the experiment proposed for flight. 
It is also an EP experiment in its own right aiming to 
compete with torsion balance tests 0,0. ^ n this effort, 
motor noise, low frequency terrain tilts |16| and tidal per- 
turbations |17| are the main issues to be addressed. 

A full knowledge of the dynamical response of the GGG 
rotor is needed, especially in view of its condition of su- 
percritical rotation, and of its common mode rejection 
behavior. The theoretical understanding of the dynami- 
cal properties of the rotor, together with the construction 
of a full simulation facility, would allow us to predict and 
interpret the collected experimental data; they also pro- 
vide a virtual environment for planning the experiment 
and optimizing its performance. 

With these motivations in mind, we demonstrate that 
a simple but very effective mathematical model can be set 
up to quantitatively describe the dynamical properties of 
the GGG rotor. In this paper (Part I), we determine the 
normal modes in all regimes, from subcritical to super- 
critical rotation, and address the issue of self-centering in 
supercritical rotation. In the following paper (Part II), 
we provide the dependence of the common mode rejec- 
tion ratio on various system parameters which govern the 
design of the instrument. 

The differential equations in the model are solved by 
means of a user-friendly simulation program and the nu- 
merical solutions are tested against the data available 
from the experiment. The physical content of the model 
is also discussed by means of approximate analytical so- 
lutions, which provide useful physical insight. 



This paper is organized as follows: Sec. |H] describes 
the main features of the experimental apparatus; Sec. 11111 
presents the dynamical model of the system, referring to 
specific appendices for details. Sec. IIVI reports on the 
numerical method that we have implemented and Sec.lVl 
gives the results obtained on the determination of the 
normal modes of the system, showing an excellent agree- 
ment between theoretical predictions and experimental 
data. The details of the calculations are contained in two 
appendices, while the third one is specifically devoted to 
the important concept of self-centering. Concluding re- 
marks are discussed in Sec. IVII 



II. THE GGG ROTOR: OVERVIEW OF THE 
EXPERIMENT 

GGG is a rotating differential accelerometer operated 
in a vacuum chamber (see Fig.^l. It is made of two con- 
centric hollow test cylinders, 10 kg each, weakly coupled 
by means of a vertical arm —a tube located along the 
axis of the cylinders— to form a vertical beam balance 
(from now on we shall always omit the term "hollow" 
when referring to the test cylinders). The coupling arm 
is suspended at its midpoint from a rotating vertical shaft 
in the shape of a tube enclosing it (see Fig. ^ and Fig. 
right hand side) . A total of three suspensions are needed 
(drawn in red in Fig. : a central one (see Fig. |2 left 
hand side) to suspend the coupling arm from the rotating 
shaft, and one for each test cylinder to suspend each of 
them from the top and bottom end of the vertical cou- 
pling arm. 

The suspensions are cardanic laminar suspensions 
manufactured in CuBe which are stiff in the axial direc- 
tion Z , against local gravity, and soft in the plane X — Y 
orthogonal to the axis so that the geometry is naturally 
two-dimensional, the horizontal plane being sensitive to 
differential accelerations acting between the test cylin- 
ders. In normal operation mode, modulation of such a 
signal is provided by setting the whole system in rotation 
around the vertical axis in supercritical regime, namely 
at frequencies v s larger than the natural differential fre- 
quencies of the rotor, typically v s > 1.5 Hz. The dif- 
ferential character of the instrument is strengthened by 
two differential read-out systems made of four capaci- 
tance plates (indicated as IP, internal plates, in Fig.^) 
located in between the test cylinders and which are part 
of two capacitance bridges in two orthogonal directions 
of the sensitive plane. 



A. Description of the mechanical structure of the 
apparatus 

The GGG apparatus is schematically presented in 
Fig.^ where a section through the spin- symmetry axis Z 
is shown inside the vacuum chamber C. At the top-center 
of the frame is the motor M whose shaft is connected to 
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FIG. 1: Section through the spin axis Z of the differential 
accelerometer inside the vacuum chamber. C: vacuum cham- 
ber; M: motor; OD: optical device (see Sec. Ill Bfl : x: ball 
bearings; ST: suspension tube; A: coupling (balance) arm, 
located inside the suspension tube, with its 3 laminar cardanic 
suspensions (in red); G: center of mass of the two cylinder's 
system (in blue the outer cylinder, in green the inner one, 10 
kg each). IP are the internal capacitance plates of the dif- 
ferential motion detector (Sec. Iff 51 . OP are the outer ones 
for whirl control (Sec. Ill C 2L PC is the contactless inductive 
power coupler providing power to the electronics inside the ro- 
tor. The relevant distances Li, and L of the centers of mass 
of the inner and outer bodies from their suspension points are 
also sketched, along with the arm length 2L a + AL. T and P, 
at the top of the rotor, are the tiltmeter and 3 PZTs (at 120° 
from one another -only one shown) for automated control of 
low frequency terrain tilts. The drawing is to scale and the 
inner diameter of the vacuum chamber is 1 m. 

the suspension tube of the rotor ST (drawn in yellow) 
by means of an appropriate motor-rotor joint and turns 
in the vertical direction inside ball bearings, indicated by 
x symbols in the figure. From the suspension tube ST 
rotation is then transmitted to a tube located inside it 
which constitutes the vertical beam of the balance (also 
referred to as the coupling arm, Fig. El right hand side), 
the connection between the two being provided at the 



midpoint of the arm by the central laminar cardanic sus- 
pension (see Fig. |21 left hand side) . 

The coupling arm in its turn transmits rotation to both 
the test cylinders, since they are suspended (by means of 
two laminar cardanic suspensions similar to the central 
one) from its two ends. More precisely, the inner test 
cylinder (shown in green) is suspended from the bottom 
of the coupling arm at a distance Li from the cylinder's 
center-of-mass, while the outer one (shown in blue), is 
suspended from the top of the coupling arm, at a dis- 
tance L Q from the cylinder's center-of-mass. In Fig. 
the three suspensions are drawn in red. It is apparent 
that the central suspension carries the whole weight of 
the rotor, mostly the weight of the two test cylinders 
(10 kg each) plus the small weight of the coupling arm. 
It is worth noting that the metallic suspensions provide 
passive electrostatic discharging of the test masses. 




FIG. 2: Left hand side: the central laminar cardanic sus- 
pension of the GGG rotor, located at the midpoint to the 
coupling arm in order to suspend it from the suspension tube 
(shaft). Right hand side: the coupling arm inside the suspen- 
sion tube (shaft) as seen from the top. Two cardanic laminar 
suspensions are located at its top and bottom ends. They sus- 
pend the test cylinders (not shown here) through two metal 
rings. The dimensions of the rings depend on the dimensions 
of the concentric cylinders, which have equal mass (10 kg) and 
therefore different size. Top and bottom ring refer to outer 
and inner test cylinder respectively. 

In this way, the symmetry of the whole apparatus is 
cylindrical, its axis being both the vertical beam of the 
balance and the axis of rotation, the balance is sensitive 
in the horizontal plane and the test masses are concentric. 



B. The differential motion detector system 

The Differential Motion Detector (DMD) reflects the 
cylindrical symmetry of the system and is composed of 
the following three parts. 

(1) Two capacitance plates IP working as DMD(x) 
sensors (drawn as vertical lines in between the cylinders 
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in Fig.^l, are located halfway in between the test cylin- 
ders in correspondence to the X direction with a clear 
gap of 5 mm on either side and are connected to the 
suspension tube by means of an insulating frame. A sim- 
ilar pair of capacitances forming a DMD(y) is placed in 
the Y direction. A voltage signal is applied to each ca- 
pacitance bridge in order to shift the signal of interest 
to a high-frequency band with reduced 1// noise (with 
phase locked detection). The filtered signal is digitized 
by an ADC before transmission to the non-rotating (lab- 
oratory) frame. Calibration and balancing of the capaci- 
tance bridge are performed by means of procedures out- 
lined in • The best sensitivity achieved in bench tests 
corresponds to mechanical displacements of 5 pm in 1 s 
of integration time 0, ^| . Presently, the sensitivity of 
the read-out system during normal operation is ~ 10~ 9 
m. 

(2) An optical device OD located below the motor and 
above the ball bearings, utilizing a disk with 32- holes and 
an infra-red emitter-detector pair, provides a reference 
signal for the angular position of the rotor. The refer- 
ence signal is combined with the X and Y channel data 
from the DMD and encoded into RS232 format for trans- 
mission to a computer. Then a second emitter-detector 
pair located at the very bottom of the rotor (using a hole 
along the axis of the power coupler PC, see Fig.^l trans- 
mits the digital signal from the rotor to the non-rotating 
frame from where it is taken out of the vacuum chamber 
through electrical feed-throughs. 

(3) An annular disk, in two semicircular parts, is 
mounted around the upper half of the suspension tube 
and contains the two capacitance-bridge circuits and 
their preamplifiers (see Fig.^l. The necessary electronics 
to demodulate the signal and convert it from analogue to 
digital form, as well as the drivers for the optical emitter 
are also located here. 



C. Principle of operation 

For detecting an EP violation signal the instrument 
relies on its sensitivity to relative displacement of the 
two test masses, which in the final design will be made 
of different materials. An acceleration in the horizontal 
plane of the laboratory acting differently between the test 
cylinders gives rise to a relative displacement of the two 
in the direction of the acceleration. This displacement 
unbalances the capacitance bridges and gives rise to an 
electric voltage proportional to it. 

A modulation of the displacement, as seen by the ca- 
pacitance plates, is achieved by setting the whole system 
in rotation around the vertical axis of symmetry passing 
through the shaft, as shown in Fig.Q] Note that the sig- 
nal modulation obtained in this way does not affect the 
centers of mass of the test cylinders, hence it does not 
affect their relative displacement, which is the physical 
quantity measured in the experiment. As a result, this 
type of modulation reduces the noise but not the signal. 



In point of fact, this signal modulation could be 
achieved by keeping the test cylinders stationary and ro- 
tating only the capacitance plates (located in between 
the two, indicated as IP in Fig.QJ which form the differ- 
ential motion detector system described above. However, 
by rotating the test cylinders together with the capaci- 
tors, any irregularity in their mass distribution averages 
out; moreover, the supercritical regime can be exploited 
to reduce rotation noise for all parts of the apparatus 
(see Sec. Ill C "2*1 and Appendix [UJ. As for the experiment 
in space, rotation of the whole spacecraft has two more 
very important advantages. In the first place, it elim- 
inates the need for motor and ball bearings altogether, 
which are a considerable source of noise in the ground 
experiment. Secondly, by rotating around the axis of 
maximum moment of inertia, the spacecraft is passively 
stabilized, thus reducing its weight, cost and complexity, 
as well as disturbances on the EP experiment. 

An EP violation signal in the gravitational field of ei- 
ther the Earth or the Sun would have a component in 
the horizontal plane of the laboratory which could be de- 
tected by the instrument. Since the test bodies are rotors 
suspended on the Earth, and the Earth rotates around its 
axis, this diurnal rotation gives rise to large gyroscopic 
effects on the test bodies resulting in a non zero differen- 
tial acceleration which would mask an EP violation signal 
in the field of the Earth itself. Measurements of such gy- 
roscopic effects have been reported in (Sec. 5, Fig. 
12). The instrument —in this ground based version— is 
therefore used for two purposes: i) to establish its sensi- 
tivity as a prototype of the flight instrument, namely for 
an expected signal at the orbital frequency of the satel- 
lite (~ 1.75 • 10~ 4 Hz, i.e. about 1 and half hour period, 
at an Earth orbiting altitude of ~ 520 km); ii) to look 
for an EP violation in the gravitational field of the Sun, 
in which case the signature of the signal (see 0, Sec. 
2) would have a dominant Fourier component of 24-hr 
period due to the diurnal rotation of the Earth. 



1. Differential character and common mode rejection 

The differential character of the whole instrument, 
namely its capability to reject accelerations which are 
common to both test masses, is in principle ensured by 
the geometry and mounting of the test masses. It is fur- 
ther augmented by the differential nature of the DMD 
system. 

(i) The sensitivity of the instrument to differential ac- 
celerations of the test masses depends on the softness of 
the laminar suspensions and on the uniform distribution 
of mass around the spin axis. Soft suspensions and a 
good balancing of the rotor provide long natural periods 
for differential oscillations of the test masses relative to 
each other, giving rise to larger relative displacements 
between the two, and in turn to stronger output voltage 
signals. 

Tuning of the natural differential period Try of the test 
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cylinders is made possible by changing a moment-arm 
in the beam balance. This is accomplished by moving a 
small solid ring mounted at the lower end of the balance 
(coupling) arm. Moving this ring vertically along the 
arm, in the Z direction, displaces the center-of-mass of 
the balance arm from its suspension point by a quantity 
AL. If AL = the center-of-mass of the balance arm 
is coincident with its suspension point. AL can be ad- 
justed to be either slightly positive or negative, resulting 
in a longer or shorter Tr>. However, there's a maximum 
positive value that AL can assume before the system be- 
comes unstable (see Eq. (|4*T]) below). 

Asymmetric distribution of mass of the rotor in the 
horizontal plane, resulting in a non-zero inclination of 
the coupling arm in the rotating reference frame, may 
also be corrected by two small masses mounted inside 
the coupling arm itself, one of which is movable in the X 
direction and the other in the Y. 

The tilt of the spin axis with respect to the non rotat- 
ing laboratory frame is controlled by three micrometer 
screws which support the plate on which the rotor shaft 
is mounted. In addition the tilt can be finely adjusted 
using piezoelectric actuators (P) attached to the tips of 
the micrometer screws (see Fig.0. 

(ii) As to the DMD system, a non-zero off-centering 
of the capacitor plates IP located in between the test 
cylinders —measured by the ratio (a — b)/a where a (b) is 
the nominal gap between the inner (outer) mass and any 
one of the capacitance plates— would make a common 
mode displacement Axe of the test masses to produce 
a differential output signal in addition to that produced 
by a real differential displacement Axr>. The larger this 
off-centering, the larger the fraction of common mode dis- 
placement which is turned into a 'fake' differential signal, 
i.e. which contributes to the total unbalance AC of the 
capacitance bridge [llj (GG Phase A Report, Sec. 2.1.3) 
from the original capacitance value C : 

AC a-b A 1 

T^r - -^ A *c - —Axd ■ 1 



2. Signal modulation and whirl motions 



Hz) higher than the natural frequencies v n of the system, 
a condition known as supercritical rotation. 

The GGG apparatus has 3 natural frequencies. The 
differential frequency vry, of the oscillations of the test 
bodies relative to one another, and two common mode 
frequencies, vq\ and vc2, of both test masses together. 
In the GGG setting reported here their values are: Vd = 
0.09 Hz, v C \ = 0.91 Rz^c2 = 1-26 Hz. 

It is well known 0, EoL l2l| that in supercritical rota- 
tion the masses are able to self-center and greatly reduce 
the original offsets of their centers of mass with respect 
to their own rotation axes. Any initial offset, which in- 
evitably results from construction and mounting errors, is 
in fact reduced by a factor {vd/^s) 2 - Such self-centering 
is a very essential requirement when using fast rotating 
macroscopic test bodies for the purpose of detecting the 
effects of extremely small forces between them. 

It is also well known that in supercritical rotation, dis- 
sipation in the system gives rise to destabilizing whirl 
motions at frequencies v w equal (or close) to the natu- 
ral frequencies of the system, whose amplitude increases 
with time at a rate 1/t w = ~nv w j 'Q(v s ) scaling as the 
whirl frequency v w and the inverse of the quality factor 
Q at the spin frequency v s [21II22I l23l| . 

Whirls can be stabilized by passive and active meth- 
ods. Passive stabilization is typically used in engineering 
applications of supercritical rotors, but it produces too 
large disturbances for our purposes. We have used a pas- 
sive damper in the past only to stabilize the rotor during 
resonance crossing (see [l4(, Sec. 3). With the current 
improved apparatus, damping at resonance crossing is no 
longer needed. A much finer whirl stabilization can be 
performed actively by means of 8 small capacitance sen- 
sors/actuators (indicated as OP, outer plates, in Fig.^) 
placed close to the outside surface of the outer test cylin- 
der, 4 of them used as sensors and 4 as actuators in two 
orthogonal directions of the horizontal plane . In the 
GGG experiment performed at supercritical speed the 
relevant Q value is determined by losses due to deforma- 
tions of the laminar suspensions at the spin frequency. 
Experimental measurements of Q are reported in [l4j 
and, more recently, in [l8| 



Signal modulation in testing the equivalence princi- 
ple has been first proposed in 2] in order to improve 
the Eotvos experiments. By referring to the Sun rather 
than the Earth as the source mass of the gravitational 
field, the diurnal rotation of the Earth itself on which 
the test masses are suspended provides a 24-hr modula- 
tion with no need to rotate the experimental apparatus, 
an operation which gives rise to relevant disturbances in 
such small force experiments. However, higher modu- 
lation frequencies are desirable in order to reduce 1// 
noise, and in fact excellent results have been obtained 
by 0, with torsion balances placed on a turntable ro- 
tating faster than the Earth. In GGG we try to spin the 
test masses much faster, at frequencies (typically a few 



III. THE MODEL 

Having described the real instrument, we are now in a 
position to outline the minimal model used to describe 
its dynamical behavior. Fig. [21 displays a schematic rep- 
resentation of the model in the reference frame {X'Y'Z'} 
rotating with the shaft at an angular velocity lo s = 2-kv s 
around the Z' axis (cD s = lo s Z'). The relevant parts of 
the instrument depicted in Fig. ^ are sketched in Fig. 
with the same colors. The coupling arm, with mass m a 
(drawn in cyan as in Fig. and length 2L a + AL is 
suspended at its midpoint MP from the rotating shaft 
and suspension tube ST (yellow) by means of the cen- 
tral laminar suspension LS (red) with elastic constant 
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K. The vector e is the offset of the arm center-of-mass 
from the axis, which is unavoidable because of construc- 
tion and mounting errors. Variations of AL, as we have 
already discussed, produce a change of the mass distribu- 
tion, hence of the natural differential period of the test 
masses, To. Here, and with no loss of generality, e is 
placed along the X' axis. 

The outer test cylinder, of mass m (blue), is sus- 
pended from the top of the coupling arm by means of 
the laminar suspension with elastic constant K a and 
its center of mass is at a distance L a from the sus- 
pension. In a similar manner, the inner test mass rrii 
(green) is suspended from the bottom of the arm, Ki 
and Li being the corresponding parameters. From now 
on, the label A = z,o, a will be used to refer to the pa- 
rameters of the inner mass, outer mass, and coupling- 
arm respectively. The three bodies have moments of 
inertia 7 A £ = ^Ar, = m x {3R 2 XI + 3R 2 XE + R 2 XH )/12 and 
I\( = m \{R\i + R\e)/^ along their principal axes (£, -q 
and C) , R\j , R\e being the internal and external radii of 
the cylinder A, and R\h its height. 




FIG. 3: Minimal model for the real instrument sketched in 
Fig. Q (see text for details) . On the left hand side the var- 
ious parts are drawn with the same colors and labels as in 
Fig. Q Here the midpoint of the coupling arm is indicated 
as MP. L a , L and Li refer to the dimensions of the cou- 
pling arm and the outer mass and inner mass suspension arms 
respectively. L a , L , and Li are the unit vectors of the cor- 
responding beams. The offset vector e, due to construction 
and mounting imperfections, is also indicated. On the right 
hand side we sketch one of the cylinders in the rotating ref- 
erence frame {X'Y' Z'}, showing its principal axes of inertia 
{Ci^jCji the position vector R\ (A = a,o,i) of its center of 
mass and the angles 8\,<f>\, which are not the usual Euler an- 
gles, as discussed in the text. Below this figure, the small one 
to the right shows a typical deformation of one of the laminar 
suspensions of length for instance the central one. None of 
these figures is to scale. 

The laminar suspensions have length I, the central 
one is slightly stiffer than the other two and we assume 
Ki = K a 7^ K. In a refined version of the model, and 
whenever specificed, we also consider an anisotropic cen- 
tral suspension by introducing the parameter A such that 
K Y '=AK X >. 



By defining the unit vector L a of the coupling arm as 
pointing from its midpoint towards the bottom suspen- 
sion, and the unit vectors L a and Li of the test cylinders 
each pointing from the suspension to the center of mass 
of the body (see Fig. O , the corresponding position vec- 
tors in the rotating reference frame {X'Y' Z'} of Fig. [3J 
are 

R a = e-0.5ALL a 

Ro = e-(L a + AL)L a + L L (2) 
Ri = e + L a L a + LiLi . 

A. The Lagrangean 

The Lagrangean C in the rotating reference frame 
{X'Y'Z'} can be written as 

C = T-U, (3) 

where the kinetic term can be very generally written as 

T =\zZj vldmx (4) 

after defining the velocity v\ of the mass element dm\ 
in body A with volume t\. Then, U includes the poten- 
tial energies associated with gravity and with the elastic 
forces, namely 

U = U g + U el (5) 

where 

U g = ^ -mxg-Rx with g=-gZ' (6) 

A— a,o,i 

U el = \ K * l2 \R* x i? A | 2 (7) 

X—o.i 

+ h 2 a (K x ,\R a X'\ 2 + K Y ,\R a Y'\ 2 ^j . (8) 

For the expression of U e i we refer to the small figure at the 
bottom right of Fig. [21 sketching the laminar suspension 
and its orientation. 

We proceed along the main steps to derive the opera- 
tional expression for C The bodies are rotating around 
their own axis with angular velocity oj s in a reference 
frame which is rotating as well, as sketched in Fig. [31 right 
hand side. Thus, we define as (l\ the angular-velocity 
vector of the element dm\ in the {X'Y 1 Z'} frame and 
lu s = lu s Z' , so that 

v x = V x + A x fx = V\ + fi A x {Rx + pa) , (9) 

where Vx is the velocity of the center-of-mass of body A 
and fx is the vector pointing to the element dmx, com- 
posed by Rx and px as drawn in Fig. [31 By inserting 
Eqs.© into Eq. |0J), we can write T as 

T = T kin + Tear + Ucor + U c , (10) 
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where the only non-zero terms are (see Appendix ^ for 
details) 



Tcor + U( 



Cor 



_ ( n ) 

X m A V? A ■ (w s x i? A ) (12) 

A 

x pa ) • (w a x p A ) dmXl3) 



A Jr A 
Z A J ^ 



dm\(14) 



In Eqs. p2|) - (|13|) the terms coming from Coriolis forces 
have been split into the Ucor potential energy, which con- 
tains only the position vectors, and Tcor which contains 
also the velocities. The centrifugal part U c has been indi- 
cated as a potential energy. To proceed further, we now 
have to specify the choice of the generalized coordinates. 



B. Choice of the generalized coordinates 

The GGG rotor model shown in Fig. |3] is composed 
of n& = 3 coupled bodies, for a total of 18 degrees of 
freedom. However, the central suspension prevents them 
from performing translational motions, thereby reducing 
the degrees of freedom to 9. In addition, the motor forces 
the three bodies to rotate at a constant angular velocity, 
so that the number of degrees of freedom for the model 
is n = 6. 

We have chosen as generalized coordinates for each 
body the two angles 9\ and <fix (see Fig. right hand 
side). These angles are defined slightly differently from 
the usual Euler angles: 9\ is the angle between R\ and 
the axis —Z' and runs in the interval [0, tt]; </> a is the 
angle from the X' axis to the projection of R\ on the 
X'Y' plane and runs in the interval [0, 2tt). We thus de- 
fine the vector Q of the generalized coordinates and the 
corresponding velocities Q 

Q = {(11,12, <?3, <74, <?5, Qe} = {0a, <Pa,do, <j>o, h, <M ■ (15) 

With these definitions in hand, we have that 

L a = {sm9 a cos0 Q , sin6> Q sin</> a , - cos6> a } , (16) 

and similar expressions for L a and Li. Eq. Ill Gil turns Eqs. 
(J2J) into expressions for the R\(Q) and the corresponding- 
velocities V\ = R\(Q,Q). We then conveniently write 
all the vectors in the {X'Y'Z 1 } reference frame in terms 
of their components in the {^C} frame by means of the 

rotation matrix M fEq. lA2| ). namely 57 A =jv[ ^A.f^c and 

P\ =M P\,£riC- 

After noting that Q\ } £ri( — {~ &\, 0a sin# A , uj s } and 
performing all the integrals over the three bodies, we fi- 
nally obtain (Appendix ^ the operative expression for 



C(Q, Q) in the rotating reference frame: 

C(Q,Q)=T(Q,Q)-U(Q), (17) 

where we have defined 

T(Q, Q) = T kin (Q, Q) + T Cor (Q, Q) (18) 
U{Q) = U g (Q) + U e i(Q) - UcoAQ) - U C (Q) . (19) 

The terms entering <|17f) - (ll9fl are 

T km = [mxV x (Q,Q) 2 +ht(<Plsm 2 8 x + 6 

A 

Tear = ^m A yx(g,Q) • (cT, x R X (Q) 



(20) 



+ X I a£W s ^a sin 2 9 x , 



U C or = - X hc^s cos 9 X , 

A 

and 

U c = -^2m x (u s x R\(Q) 

A 



\ ^ [hi sin 2 9x + he cos 2 Ox] u 2 s . (21) 



To these equations we have to add the expressions ©-0 
written in terms of Rx{Q) through (J2J and (| 1 . 

Eqs. together with Eqs. (19)-J2TJ yield the La- 
grange function of the model in Fig. [21 



C. Equilibrium positions and second-order 
expansion 

During normal and successful operation of the GGG 
rotor only very small amplitude motions take place. The 
Lagrange function (|17f) can thus be expanded to sec- 
ond order in (Q,Q) around the equilibrium solution 

(Q°,Q — 0), Q° = {q®, ...,96} to derive linearized equa- 
tions of motion. 

In order to do this, we first determine the equilibrium 
positions from the equation 



dU 

dq~^= q ° 







1. 



(22) 



We then use the physical assumption that during the 
motion, the Q's are sligthly perturbed from their equi- 
librium values Q . This results in the substitutions 



Q 
Q 



Q° + Q 
Q. 



(23) 
(24) 



into l|17H to obtain a linearized version of the Lagrange 
function. £(Q, Q) can now be expanded to second order, 
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namely 

n n n 

C(Q,Q) = a + ^a jk q J q k + ^b jk qjqk + ^ c jk qjqk 

j<k j<k j,fc=l 



^djqj +0(q jl q k ) 4 



(25) 



where we remark that now the g^'s are small according 
to the substitutions 1(23 - 1(24 |) . and that the linear terms 
have cancelled out because of (12211 . The matrix coeffi- 
cients djk, bjk, and Cjk are known functions of the Q° 
and of the governing parameters of the system, and in 
general are to be numerically evaluated. 



D. Linearized equations of motion 

The equations of motion in terms of the known djk, 
bjk, and Cj k coefficients are 



d dC dC 



j = 1, ...,n 



dt dqj dqj 
where we have introduced the generalized forces 

3 



?i - ^2 Fxd 



dRxd 
d qj 



(26) 



(27) 



starting from the cartesian components F\d of the forces 
acting on each body. The Tj are to be consistently ex- 
panded to first order, namely 



Tj = ^a jk qk+^(3 jk q 



Qk 



(28) 



k=l 



k=l 



By combining Eqs. 1)25 (1 - l(2ij |) together, the equations of 
motion can be written in a compact matrix form as 



= S 



(29) 



with the obvious notation Q = {q\, ...,qe}. In Eq. 1(29(1 . 

M is the n x n (n = 6) "mass"-matrix composed by the 
bjk coefficients: 



l jk = 2b lk 6 lk + b jk (l-Sjk) , 



(30) 



where the factor of 2 on the diagonal elements is a conse- 
quence of the restricted j < k sum in the expansion (|25|l . 
§ is a n x 2n matrix containing the djk, Cjk, otjki an d Pju 
coefficients: 



§ = A 2 + C 2 + Ai + Bi , 



(31) 



with 



\2jk=^2k] = Za^Sjk + a jk (l - 5jk) k<n (32) 
= n < k < 2n , 



C 



2jk 



^2kj — Cjk — Ckj 



\l k j — Ujk 

= 



and 



Hjk 



nkj 



Pjk 



k<n (33) 
n < k < 2n , 

k < n (34) 
n < k < 2n , 



k < n (35) 
n < k < 2n . 



Note that while M and the submatrix defined by the first 
ri = 6 columns of A 2 is symmetric, the submatrix defined 
by the second n = 6 columns of C 2 is antisymmetric, as 
expected after inspection of the expansion (|25|l . 

For all practical purposes, it is convenient to turn 129(1 
into a more symmetric form involving only first-order 
time derivatives. To this aim, we define the 2n = 12- 
component vector X as 



X2j-l - <lj 



j = l, ...,n = 6 (36) 



X 



2j - 1j 



By inserting the definition 1(36(1 into 1(29(1. we finally ob- 
tain 



X = AX , 



(37) 



where A is now the square 2n x 2n dynamical matrix 
defined from M _1 and § after inserting rows of zeros. 



.4 



(M- X S) 



j even and k odd 

j even and k even 

j odd and fc = j + 1 
j odd and fc ^ j + 1 



(38) 



The relations ((37(1 - 1(38(1 are central equations, written 
in a form amenable for numerical evaluation. The eigen- 
values of the A matrix 1(38(1 correspond to the normal 
modes of the rotor, and the solution of the set of dif- 
ferential equations ((37(1 completely determines the small- 
amplitude dynamical behavior of the rotor modelled in 
Fig- El Before turning to the description of the numerical 
method, we introduce rotating and non-rotating damp- 
ing. 



1. Rotating and non-rotating damping 

By means of (I27H - ((28(I we can in principle introduce 
any known force determining the dynamical behavior of 
the rotor. In the following we include dissipative forces 
TZr and IZnr due to rotating and non-rotating damping 
mechanisms respectively (see [2(J,|2l|]). The rotating part 
of the dissipative force is to be ascribed to dissipation of 
the laminar suspensions. In supercritical rotation, this 
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kind of dissipation is known to destabilize the system, 
generating whirl motions. It can be expressed as 

n R {Q,Q) = -T Ra (L a - (L a ■ Z') ■ Z') (39) 

- Y, r ^(^ - • Ra(Q))Ra{Q)) , 

A — 0,2 

where the velocities L\ are functions of {Q,Q). Instead, 
non-rotating damping has the effect of stabilizing a sys- 
tem in supercritical rotation, and can be written as 

n NR (Q,Q) = -T NR (Vi(Q,Q)+ti s x R^Q)) . (40) 

Other forces acting on the rotor, such as external dis- 
turbances due -for instance- to tides and seismic noise, 
or control forces applied in order to control the rotor dy- 
namics, can also be included, as described in Sec. IIVI 

IV. THE NUMERICAL METHOD 
A. General considerations 

The simulation method that we have implemented, rig- 
orously follows the derivation outlined in Sec. 11111 We 
have found very convenient to use the MATLAB envi- 
ronment, with SYMBOLIC TOOLBOX and SIMULINK 
packages, as it allows us to perform all the needed sym- 
bolic calculations and numerical evaluations, together 
with the analysis of experimental data. 

We start from the formal Lagrange function written in 
a user- friendly way as in |JS}-|JSJ and i|llll4[l by means 
of symbolic vector operations. We specify the choice 
(|15ll for the generalized coordinates with respect to the 
{X'Y'Z'} reference frame and define accordingly all the 
vectors entering C We then move on to the symbolic 
computation by linearizing and expanding the Lagrange 
function as in l|25l) . and define the matrices M, A2, C2, 
Ai, Mi, and A. 

Once the system parameters are fixed (see below), 
the numerical computation is carried out using standard 
packages to find eigenvalues and eigenvectors of the A 
matrix, that are the normal frequencies and modes of 
the spinning rotor. The A matrix is then inserted as 
input to perform the dynamical simulation within stan- 
dard transfer-function method used in the SIMULINK 
toolbox. 

The advantage of this strategy is apparent, in that it 
easily allows us to make any changes in the model that 
correspond to changes in the experiment we would like to 
test before implementation. Since the number of bodies 
rib and of the generalized coordinates n are symbolically 
defined and specified only once, all what is to be done in 
order to introduce any changes or new features amounts 
to modifying or add pieces of the Lagrangean after having 
symbolically written them in terms of vector operations. 

The description of the method used to introduce exter- 
nal forces is postponed to Part II of this work, where it 



is used to evaluate the common mode rejection function. 
We now turn to listing the system parameters. 



B. System parameters 

The parameters which govern the physics of the GGG 
rotor are: the geometrical dimensions of the three bodies, 
their weight, the mounting error e, the elastic constants, 
length and anisotropy factor A of the three laminar sus- 
pensions, the quality factor Q. To these parameters - 
—which are fixed after construction— we must add the 
spin frequency v s — lo s /2tt, that can be varied in the 
course of the experiment. The balancing of the beams 
and the natural period Tjj of oscillation of the test cylin- 
ders relative to one another can also be adjusted as dis- 
cussed earlier by moving small masses along the balance 
(coupling) arm. 

We have inserted as inputs to the numerical calculation 
all the above parameters as determined in the real GGG 
instrument. They are listed in Tab. HI and HJ As for the 
spin frequency, in the experiment it varies in the range 
(0 < v s < 3.9) Hz while in the model calculations it can 
be assumed in a wider range (0 < v s < 10) Hz. 

The differential period Tr> corresponding to the value 
AL listed in Tab-HJis measured to be 11.7 s and 10.8 s in 
the X and Y directions respectively. These values are in 
reasonable agreement with the following simple formula: 



/ (K+Kj+K„)l 2 g AL 
V (m;+m )L2 2L a L a 

which can be derived from the general equations of mo- 
tion (|26|l describing the small oscillations of the 0\ angles, 
in the very simplified case in which the bodies are nei- 
ther rotating nor subjected to any dissipative or other 
external forces, except gravity, and under the reasonable 
assumption that Qi = o = and that <f>\ 's are costant, 
e.g. (j>\ = 0. 

Of the whole set of parameters used, only the 
anisotropy factor A of the suspensions and the construc- 
tion and mounting error e are not measured from the 
instrument. A is tuned, together with the balancing AL, 
so as to reproduce the natural frequencies of the non- 
spinning instrument. A conservative value e = 20 /ini is 
assumed for the offset, and it is checked a posteriori not 
to have any sizable effect on these results. 

V. RESULTS: THE NORMAL MODES 

We have solved for the eigenvalues v n of the matrix A 
in 1(3*71) . using the system parameters listed in Sec. HVBl 

Fig. [J] summarizes our results by plotting the normal 
modes of the system (y n in the non-rotating frame) as 
functions of the spin frequency v s = lu s /2tt of the ro- 
tor. In this figure, theoretical results for v n (v s ys are 
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TABLE I: Input parameters for the numerical calculations: 
geometrical dimensions of the real bodies. (A mounting error 
of e = 20 )im has also been used) 
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Body 



Rx 



Arm (a) 



(kg) (cm) 
0.3 3.3 



(cm) (cm) 
3.5 19 

(AL = -0.1353) 
Outer cyl. (o) 10 12.1 13.1 2L a + Lt + AL 29.8 
10.9 4.5 



R\H 

(cm) 

2L a + AL 



Inner cyl. (i) 10 8.0 



21.0 



TABLE II: Input parameters for the numerical calculations: 
laminar suspensions data. In addition, a conservative value 
Q{v a ) — 510 has been used taken from previous measurements 
of whirl growth ( 14], Fig. 7). 

Suspension I K Anisotropy 

(cm) (dyne/cm) (A = K Y >/K X ') 

Central 0.5 10 a {l/L a ) 2 2.6 

Outer cyl. (o) 0.5 10 6 1.0 

Inner cyl. (i) 0.5 10 6 1.0 



displayed by the solid lines in the case of zero rotating 
damping (i.e., no dissipation in the suspensions), and by 
the open circles in the case of non-zero rotating damp- 
ing (Q(v s ) = 510). At zero damping there are 12 lines, 
6 horizontal and 6 inclined; starting from the 3 natural 
frequencies (Sec. Ill 0"2)l of the system, we get 3x2x2=12 
normal mode lines, a factor 2 being due to anisotropy 
of the suspensions in the two orthogonal direction of the 
plane, and the other to the positive and negative sign 
(i.e. counter-clockwise or clockwise whirl motion). For 
the non-zero damping case (open circles) a conservative 
low value Q(v s ) — 510 has been assumed (see Tab. ITT)l . 
referring to a comparatively large dissipation. This value 
has been obtained from previous not so favorable mea- 
surements of whirl growth, while much higher Q values 
(namely, much smaller dissipations) are expected (see 
discussion on this issue in [13, Sec 3). In Fig. H to be 
compared with the above theoretical results, we plot, as 
filled circles, the experimental results too, finding an ex- 
cellent agreement between theory and experiments. Since 
Fig.QJcontains the crucial results of this work, it is worth 
discussing its main features in detail. The main features 
are: the comparison with the experiment, the role of 
damping, the behavior at low spin frequencies, the so- 
called scissors's shape, the splitting of the normal modes 
and the presence of three instability regions. 



A. Comparison with the experiment 

In the experiment, the rotor is first accelerated to spin 
at a given frequency v s . Then, the natural modes are 
excited by means of capacitance actuators (indicated as 
OP in Fig. n in the X' or Y 1 directions at frequencies 
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FIG. 4: Normal modes of the GGG rotor: the frequencies 
of the normal modes are plotted as functions of the spin fre- 
quency u a . The normal modes as predicted theoretically as- 
suming anisotropic suspensions are shown as 12 solid lines in 
the case of zero rotating damping, and as open circles in the 
case of non-zero damping (see text). The experimental results 
are plotted as filled circles, and clearly agree with the theo- 
retical predictions. The bisecting dot-dashed line v — v s sep- 
arates the supercritical (u B > v) from the subcritical (v s < v) 
region. Three vertical thick lines are plotted in correspon- 
dence of 3 instability regions, their thickness referring to the 
width of the regions (see Sec. IV Gl 



close to the natural frequencies v° = v n (v s — 0) of the 
system at zero spin. The excitation is performed for sev- 
eral (typically ten) fundamental cycles l/i^ by means of 
voltages applied to 4 of the 8 outer plates (OP in Fig.^l. 
The actuators are then switched off, and the bodies' dis- 
placements are recorded as functions of time by means of 
the read-out described in Sec. Ill Bl Standard data analy- 
sis is then performed, by fitting the measurement data to 
extract oscillation frequencies and damping of the modes. 

Experimental data, resulting from averaging over sev- 
eral measurements, are represented as filled circles in 
Fig- El The agreement between theory and experiments is 
excellent, thus validating the model developed in Sec. lIIII 

In the experimental spectra as well as in the theoryit 
is found that the amplitudes of the modes in the subcriti- 
cal region v s < v —represented by the inclined lines, with 
their open and filled circles— are quite small, while non- 
dispersive modes (the horizontal lines, not varying with 
the spin frequency) are preferably excited. When the hor- 
izontal lines cross the inclined ones, the latter modes can 
also be excited. Since the excited modes must obviously 
be avoided in operating the experiment, this information 
is very useful, in that it is telling us that we should avoid 
to spin the system at frequencies where these line cross- 
ings occur. Even more so, spin frequencies lying in the 
instability regions must be avoided (see Sec. IV G|) 
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B. Role of damping 

We have numerically checked that dissipation present 
in the system does not significantly shift the natural- 
mode frequencies. This is apparent in Fig. ^ where the 
results obtained with damping (open circles) stay on the 
solid lines obtained in absence of damping. It is worth 
stressing the this result is especially good because, as 
discussed above, we have used a low value Q(v s ) = 510, 
corresponding to comparatively large dissipations As ex- 
pected, dissipation affects the line-shape of the peaks, 
making them wider than in absence of damping. 




2 4 6 8 10 12 5 10 15 
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FIG. 5: Normal modes of the GGG rotor. On the left hand 
side we show a zoom from Fig.^Jin the very low spin frequency 
region, showing in particular the 3 natural frequencies of the 
system in the zero spin case. On the right hand side, we plot 
a zoom from Fig.|l]in the small frequency region of both axes, 
showing the splitting into two lines of the low frequency mode 
because of anisotropy of the suspensions (the dashed line is 
the v s = v line as in Fig. |IJ. 



C. Low- frequency limit 

At zero spin frequency we have recovered theoretical 
and experimental results previously obtained for the non- 
rotating system. On the left hand side of Fig. a zoom 
from Fig. ^ at very low spin frequencies, we can see that 
the non-spinning rotor is characterized by three natu- 
ral frequencies for the instrument with ideally isotropic 
springs, the three bodies oscillating in a vertical plane. 
The frequency v a = 0.09 Hz corresponds to the differ- 
ential mode, where the centers of mass of the two test 
bodies oscillate in opposition of phase; the frequencies 
v° = 0.91 Hz and v° = 1.26 Hz correspond to common 
modes, in which the common center of mass of the two 
test bodies is displaced from the vertical. During rota- 
tion, the number of degrees of freedom increases to six, 
as discussed in Sec. IIII Bl leading to the six lines plotted 
in the same figure. 

We may get a flavor of the v s dependence of the modes 
in the v s <C v n limit, by evaluating the natural frequen- 
cies of only one spinning cylinder with mass m, moments 
of inertia 1^ and Jf. The cylinder is suspended at dis- 



tance L and with offset e from a fixed frame by means 
of a cardanic suspension with isotropic elastic constant 
K and length I. The calculation is performed by follow- 
ing the steps outlined in Sec. IIIII with rib = 1 and thus 
n = 4. After evaluating the A-matrix ( Appendix iBl see 
Eq. H-BIOI) ) and solving det(A — si) — 0, we obtain the 
two double solutions 

^1,2,3,4 = ±^ (42) 

for the four v n (y s jv n 0) in the non-rotating frame. 
In Eq. g2J>, i>° = i^y 1 ,/g/L + Kl 2 /mL 2 is the natu- 
ral frequency for the non-spinning point-like mass. L' = 
^JL 2 + (l£ - 2I c )/m (with I c < 0.5mL 2 + J 5 , see Ap- 
pendix [BJ), takes into account the extended nature of the 
body, and the ratio L/L' modifies v n with respect to v®. 



D. Anisotropy 

If the suspensions are not isotropic in the two orthog- 
onal directions, as it is indeed the case for our real car- 
danic suspensions, each natural frequency is expected to 
split up. This is clearly shown on the right hand side of 
Fig. El It is worth noting that the splitting is larger for 
the lowest frequency mode. 



E. Scissors's shape 

Fig. 21 shows that each natural frequency of the non- 
spinning system splits up into two branches at v s > 0, a 
lower branch remaining approximately constant and an 
upper branch increasing with 2v s . 

This characteristic scissors's shape can be traced back 
to general properties of spinning bodies (see also Part 
II). We again use the one-cylinder simple case (see Ap- 
pendix [0) to prove this statement. By following the 
same procedure which has led to Eq. I|42ll we obtain (see 
Eq. CE1) 



"1,2 - ±^|-2^ s -. ± , s M-^-j (43) 
^ 3 ,4 = ±\Jv 2 + 2v n v a ^ ~ ±v s [\ + V ~^\ ,(44) 

showing that in the rotating frame v n oc v s to zeroth or- 
der. After taking the v s lv n — > oo limit and tranforming 
back to the non-rotating frame by means of the substi- 
tution s n — 2mv n — + 2~K(iv n + iv s ), we finally have 

L . . 

1/2,4 =s 2v s • (46) 
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F. Mode splitting 

The two branches may cross at selected frequencies. 
Crossing and anticrossing of degenerate modes is a very 
general concept, which applies to a variety of physical 
systems, from classical to quantum mechanics, from sin- 
gle to many-particle physics. As it is well known [24| . 
splitting of the modes is expected in correspondence of 
such crossings. In our numerical results we have found 
all the 15 splittings expected for our system (see Fig.|3}. 
Fig. El shows a particular case of anticrossing of two 
modes. 




1.22- 
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FIG. 6: Normal modes of the GGG rotor. A zoom from Fig. 21 
showing one particular case of anticrossing of two modes. 



all instability regions can now be crossed, if the cross- 
ing is sufficiently fast, without producing any relevant 
disturbances even in absence of a passive damper. The 
physical space previously occupied in the vacuum cham- 
ber by the passive damper is now used for the inductive 
power coupler, indicated as PC in Fig.^ which provides 
the necessary power to the rotating electronics and has 
allowed us to avoid noisy sliding contacts. 
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FIG. 7: Normal modes of the GGG rotor. The lowest- 
frequency instability region is zoomed in from Fig. 01 



VI. CONCLUDING REMARKS 



G. Instability regions 

Dynamical instability may occur whenever the values 
of the natural frequencies are in proximity of the spin fre- 
quency. In such regions the oscillation amplitude grows 
exponentially. 

This is a well-known characteristic of rotating ma- 
chines; in engineering books it is usually described within 
the simple model of the so-called Jeffcott rotor 0] . The 
number of instability regions can be predicted from Fig. 01 
after drawing the dotted-dashed line v = v s . We have 
found indeed three instability regions. Fig. \7\ displays in 
detail the one at the lowest frequency; as shown in Fig.0] 
the two at higher frequencies are found to be wider and 
closer to each other. These theoretical results do explain 
why in the experiment we can increase the spin frequency 
and cross the low frequency instability region easily, while 
it is much more difficult to cross the frequency range 
0.9 — > 1.3 Hz. In the past we solved this problem by 
designing and installing passive dampers to be switched 
on from remote just before resonance crossing, and then 
turned off at higher spin frequencies; the least noisy was 
a special, no oil damper described in |25| . p. 45. Later on 
the GGG rotor imperfections have been reduced so that 



We have demonstrated that the linearized model set 
up in Sec. IIIII can quantitatively account for the dynami- 
cal response of the GGG rotor, an apparatus designed to 
test the equivalence principle with fast rotating, weakly 
coupled, macroscopic, concentric cylinders (Sec.^). The 
model developed here can be expanded to include exter- 
nal disturbances whose effects need to be taken into ac- 
count in testing the equivalence principle. A qualitative 
understanding has been provided, by means of helpful 
analytical solutions of the simplified model under special 
limits, of relevant features observed in the simulations as 
well as in the experimental data. 

We have acquired a detailed knowledge of the instru- 
ment's features and the way it works, the main feature 
being the normal modes of the system (Sec. [VJ in the 
whole range of spin frequencies, from subcritical to su- 
percritical regime, and as functions of the governing pa- 
rameters (see Sec lIVB)) . 

In particular: we have established the location and 
characteristics of the instability regions; we have verified 
quantitatively the effects of dissipation in the system, 
showing that losses can be dealt with and are not a mat- 
ter of concern for the experiment; we have established 
the split up of the normal modes into two scissors 's- like 
branches, distinguishing modes which are preferentially 



13 



excited (the horizontal lines) from those whose spectral 
amplitudes are typically small (the inclined lines), thus 
learning how to avoid the spin frequencies correspond- 
ing to their crossings, in order not to excite the quiet 
modes too by exchange of energy; we have investigated 
the self-centering characteristic of the GGG rotor when 
in supercritical rotation regime, gaining insight on how 
to exploit this very important physical property for im- 
proving the quality of the rotor, hence its sensitivity as 
differential accelerometer. 

In the following Part II of this work we apply the same 
model and methods developed here to investigate the 
common mode rejection behavior of the GGG rotor, a 
crucial feature of this instrument devoted to detect ex- 
tremely small differential effects. 
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APPENDIX A: THE LAGRANGE FUNCTION IN 
THE ROTATING REFERENCE FRAME 



In the following, in order to simplify the notation we 
drop the A indices everywhere, and restrict our reasoning 
to only one body. Let us begin with the expression 0}. 
After using Eq. J5} into |Q}, we have 



suits for the integrals appearing in Eq. IjAlfl : 



-2 1 - 

V dm = —mV 



dm 



-^Ti n 2 

a 

V ■ (o, x pj dm 
V ■ [ lu s x R ) dm 



V ■ (uj s x p) dm 
(l x pj ■ (l3 s x Rj dm 
x pj ■ (lu s x p) dm 
- J (lu s x Rj dm 
i / (u) s x pf dm 



~L ( 6 2 sin 2 i 



9 2 ) 

(V x Qj ■ J pdm = 
mV ■ ( uj s x R 1 



V x uj s j ■ J pdm = 
l3 s x Rj ■ J (Q x pj = 

I^u) s (j) sin 2 9 + Iquj 2 cos 9 

-m (uj s x Rj 



- [!{ sin 2 9 + I c cos 2 9] 



2 o\ , ,2 



(uj s x Rj ■ (uj s x p) dm — (ti s x Rj ■ j (uj s x p) dm = . 



By collecting all these results, one ends up with the final 
forms (19)- lPJ for the original T function 



T 



V 2 + (ti x p) 2 + 2V ■ (n x pj 



dm . (Al) 



and 



APPENDIX B: THE ONE-CYLINDER SOLUTION 



We conveniently represent the vectors fl =M. jjC 

P =M. PtvC m the HvC} frame by means of the rotation 
matrix 



(A2) 



with respect to the {X'Y'Z'} reference system. 

Thus, by exploiting the properties of the vectorial 
product and the definition of center-of-mass, namely 
J Pi = J . Pn — J Pi = 0, we find the following re- 




It is useful to study (along the lines of Sec. IIII|) the 
simplified case of only one spinning cylinder with mass 
m, moments of inertia I>t and Iq. This amounts to setting 
nt, = 1 and thus n = 4 for the number of generalized 
coordinates. 

The A matrix turns out to be 



/0 



Co 





£22 

1 





I Q.31 £23—-^14 ^33 7^-34 j 
\ £44 £44 £44 ^44 / 



(Bl) 



where the coefficients of C and 7Z are defined in terms of 
the system parameters and of the equilibrium positions 
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y , <Po as 



£11 = mui 2 L [L cos 2#o + e sin #0 ) 

- mgL cos 8 - Kl 2 cos2# + (ij - I^wf 003 26*0 

(B2) 



£-22 

£33 
£44 
£44 

•^23 

and 



I^u a cos 9q 
mL 2 + 1% 
mui 2 eL sin 6*0 
(mL 2 + / 5 )sin 2 6» 
£41 = -^c^s sin 26*o — muj s eLcosOo 
mui 2 L 2 sin 20 o 
£32 = —mu s eL cosOq , 



ftl2 
^13 
^31 
^34 



-(r fl , + rAr fl )L 2 

-rAr fl w s Lecos 6<o 
-Tnrl) s L(L sin 200 + e) 

-(r^ + rAr^^sin 2 e . 



(B3) 



(B4) 



(B5) 



In the case of negligible dissipation, the eigenvalue 
equation det(A — si) for s = 2itiv reads 

s 4 - s 2 (a 2 i + a 43 + 024042) + a 2 ia 4 3 = . (B6) 

In the uj s 3> <jJ n limit, the equilibrium solutions are 



70 



1? 



= , (B7) 



of the rotor in supercritical rotation. Let us analyze the 
one-cylinder case, by numerically integrating the equa- 
tions of motion in the presence of non-rotating damping, 
to make the rotor asymptotically stable (Sec. IIIID 1(1 . 
Fig. [3] shows the resulting motion of the cylinder in 
the horizontal plane of the rotating reference frame: its 
center-of-mass spirals inward towards an equilibrium po- 
sition much closer to the origin, i.e. to the rotation axis. 
The equilibrium position always lies in the same direc- 
tion as the initial offset vector e, which in this sim- 
ulation was assumed to be in the X' direction. The 
center of mass of the cylinder will eventually perform 
small-amplitude oscillations around the asymptotic value 
{^ = e-Lsin0 o ,r o ' = O}. 



Rotating Frame 
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FIG. 8: Self-centering of one cylinder in the presence of non- 
rotating damping. Simulated X'(t) — Y'(t) plot showing the 
motion of the center of mass of the cylinder in the rotating 
reference frame (one cylinder model, v a — 5 Hz). The center 
of mass spirals inward from the initial offset value and large 
initial oscillations to a final value, much closer to the rotation 



where V = L 2 + (7 ? - 2/ c )/m, 7 C < 0.5mL 2 + J ?) and 



a = Jl ■ KP 



(B8) 



L mL 2 



is the natural frequency of the point-like mass. Eq. I|B6|) 
becomes then 



s 4 + 2s V 



, i^ 2 



1-2 I -a 

UJ S 



(B9) 

where we have defined the natural frequency of the cylin- 
ders mass as = (LjL') 

In the LOs Wn limit, the eigenvalues equation be- 
comes instead 



L 



L 



2^(^-1 c^+l^l ^ 4 = 



(BIO) 



Eqs. ifB9)) - (jB10j) are used to derive the results 
in the main text. 



APPENDIX C: THE SELF-CENTERING 

This Appendix is devoted to a key feature of the 
GGG experiment, namely the concept of self-centering 



In the limit of small angles we obtain 

1 



4 



with 



(L'/L) 2 - (to n /io s ) 2 



</>o - (tt) 



(CI) 



(C2) 



in the case of the lower (upper) sign in Ea.f lCl|) . respec- 
tively (angles defined as in Fig.[3J). The cylinder's center 
of mass is eventually located at a distance 



AX 



±L9 () = e- 



(L'/L) 2 - {Un/LOsf 



(C3) 



from the rotation axis. 

In Fig. 13 we plot, as function of the spin frequency v s , 
the self-centering distance AX in the one-cylinder case 
discussed above, and in the point masss case. Accord- 
ing to the previous Appendix, if L' 7^ L we have the 
cylinder, while if L = L' we have the point mass. The 
two curves are worth comparing. They have a similar 
behavior till the resonance peak (in this case, at about 
v s ps 1 Hz); the distance from the rotation axis remains 
constant till, at spin frequencies slightly below the nat- 
ural one, it starts increasing showing a typical peak at 
the resonance. For the cylinder and the point mass the 
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peaks are slightly shifted. The constant value can be ob- 
tained from Eq. (|G3I) in the limit of small spin frequencies 
uo s <C u n , finding 



AX 



(C4) 



We can also recover the position and relative shift of 
the resonance peaks in the two cases from the values v v s 
(the spin frequency at the peak) taken by the poles of 
AX in Eq. (JSHJ, namely 



(C5) 



Thus, the position of the peaks is dictated by the natu- 
ral frequency v n , while the shift is due to the difference 
between L and L' . 

At rotation speeds above the resonance, and in highly 
supercritical regime u> s > w„, the behavior of the cylin- 
der and that of the point mass are remarkably different. 
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FIG. 9: Self-centering of one cylinder in the presence of non- 
rotating damping. The distance AX of the center of mass of 
the cylinder from the rotation axis is plotted as function of 
the spin frequency v s , in agreement with Eq. flCJ3ft . The same 
distance in the case of a point mass is plotted as dashed line. 
See text for comments on their comparison. 

For the cylinder, AX drops to a minimum and then 
saturates at a constant value, while for the point mass 
it keeps decreasing monotonically. The minimum for the 
cylinder is related to the presence of a zero in Eq. I|C3|I . 
namely 



±- 



L 2 



(C6) 



which is valid only if L > L. Instead, for the point mass 
we have the finite value AX = e/(l — lo 2 /lo^). Note that 
the position of the minimum shifts towards higher spin 
frequencies as L' — > L, namely, as the finite cylinder case 



approaches a point mass. In the limit oj s 3> w„ 
highly supercritical speeds) Eq. (|CJ3|> yields 



AX 



(i.e. 



(C7) 



which explains the saturation to a constant self-centring 
value in the case of a finite cylinder, whereas a point 
mass would monotonically approach perfect centering 
(i.e., AX = 0). 

In point of fact, it is very interesting to note that AX 
depends slightly on the point that we are considering 
along the cylinder's axis. In particular, in the limit w s — > 
oo, where AX of the cylinder's center of mass saturates, 
the point at distance L = L' 2 / L from the suspension 
point along the axis has instead perfect self-centering, 
namely AX = 0. This is easily seen from Eq. I|G3[) after 
substituting e ± L8q with e ± L8q and imposing AX = 0. 
We plan to exploit this property in order to obtain better 
self-centering, though it needs further investigation in the 
actual GGG rotor. 

We can consider a plot similar to that of Fig.[5|in the 
GGG case with two concentric cylinders and a coupling 
arm, where there are 3 natural frequencies (one differen- 
tial and two common mode). 

It happens that the common mode behavior is sim- 
ilar to that of the one-cylinder case (shown as a solid 
line in Fig.^J; namely, for each common mode frequency 
there is a resonance peak and a minimum peak. Instead, 
the differential frequency behavior is similar to that of 
a point mass (shown as a dashed line in Fig. [§J| . This 
latter fact is because in the differential mode the cou- 
pling arm oscillates and the cylinders' centers of mass 
move in the horizontal plane with opposite phase, while 
Oi.o = 0; under these conditions, their moment of inertia 
is irrelevant in determining the dynamics, which there- 
fore is very much alike the case of a point mass. As a 
result, the AX of the GGG rotor for intermediate val- 
ues of the spin frequency is characterized by: one peak 
at low frequency, in correspondence to the differential 
mode, and two peaks and two minima in correspondence 
to the common modes. Instead, in the limit of very low 
and very high spin frequencies, it has a behavior simi- 
lar to that displayed in Fig. depending on the values 
and directions assumed for the initial offsets of the three 
bodies. 

Thus, in order to obtain the best possible centering of 
the test cylinders in the GGG rotor, one can either spin 
at a frequency close to the minima of the common modes, 
or above both of them, in such a condition that the two 
cylinders are better centered on their own rotation axes 
than both of them are, together, in common mode. Self- 
centering on the rotation axes is very important in order 
to reduce rotation noise, because we are dealing with 
rapidly spinning macroscopic bodies and aiming at mea- 
suring extremely small effects. The issue therefore needs 
careful investigation, and to this end realistic numerical 
simulations of the apparatus are an essential tool. 
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Finally, concerning the use of supercritical rotors for 
EP testing, it is worth mentioning a frequently asked 
question: Would a relative displacement of the test bod- 
ies caused by an external force —such as that resulting 
from an EP violation— be reduced by self-centering in 
supercritical rotation as it happens for the original offset 
e? The answer is "No" , because the offset vector is fixed 



in the rotating frame of the system, while an external 
force gives rise to a displacement of the equilibrium po- 
sition of the bodies in the non-rotating reference frame. 
In the presence of such a force, whirl motion will take 
place around the displaced position of equilibrium. A 
numerical simulation, showing this important feature is 
reported and discussed in [Tl|, PLA paper, p. 176. 
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